Ring-Shaped Andreev Billiards in Quantizing Magnetic Fields 
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We present a detailed semiclassical study of a clean disk-shaped insulator-normal-metal- 
superconductor hybrid system in a magnetic field. It is based on an exact secular equation that we 
derived within the microscopic Bogoliubov-de Gennes (BdG) formalism. Results obtained from a 
classification of electron and hole orbits are in excellent agreement with those from an exact numer- 
ical diagonalization of the BdG equation. Our analysis opens up new possibilities for determining 
thermodynamic properties of mesoscopic hybrid systems. 
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Mesoscopic hybrid systems consisting of normal met- 
als (N) in contact with superconductors (S) exhibit in- 
teresting and sometimes counterintuitive equilibrium and 
transport properties resulting from the interplay between 
quantum-mechanical phase coherence and superconduct- 
ing correlations [lj, |2|. A prominent example is the 
paramagnetic re-entrance effect observed recently in ex- 
periments performed by Visani et al. on cylindri- 
cal S-Nproximity samples. While numerous theoretical 
works HJ El have addressed this problem, a fully 
satisfactory explanation of the origin of the sizable para- 
magnetic contribution to the susceptibility is still lacking. 
These previous works have studied the spectrum of An- 
dreev bound states formed in planar normal-metal layers 
in contact with a bulk superconductor, neglecting the ef- 
fects of cyclotron motion of electrons and holes due to the 
external magnetic field. Our work presented here extends 
these studies, taking into account the experimentally rel- 
evant circular geometry and fully accounting for quan- 
tum effects due to the applied magnetic field. Solving 
exactly the microscopic Bogoliubov-de Gennes equation 
(BdG)[8j, the Andreev levels in a cylindrical NS system 
are obtained for arbitrary magnetic field. In addition, we 
give a complete semiclassical description of the spectrum 
by identifying the possible classical orbits correspond- 
ing to the quantum states. This analysis is based on 
methods developed in our previous work J9| which have 
been adapted to the case of a finite magnetic field fol- 
lowing Ref. Besides being useful for shedding further 
light on causes for the above-mentioned paramagnetic 
re-entrance effect, our results are intended to serve as a 
stimulating guide to the investigation of proximity effects 
in large magnetic fields, which was the focus of several 
recent experimental an d theoretical ^3] works. 



We consider a superconducting disk of radius Rs sur- 
rounded by a normal metal region of radius Rn- (This 
models the experimentally realized 0] cylindrical geom- 
etry because motion along the axis of the cylinder adds 
only a trivial kinetic-energy term.) The magnetic field 



is perpendicular to the plane of the disk with a constant 
value of B in the N region and zero inside the S region. 
Thus, the non-zero component of the vector potential in 
pol ar coordinates (r, i?) with symmetric gauge is given by 
El A*(r,#) = B(r 2 -R 2 s )0(r-R s )/(2r 2 ), where 8 (x) 
is the Heaviside function. Excitations in an NS system 
are described by the BdG equation: 
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where VP is a two-component wave function, and Hq = 
(p — eA) 2 /(2m ff) + V — Ep. Fermi energies and ef- 
fective masses in the S and N region are denoted by 



R 



(S) T7l(N) 



E F "> and m e fi = mg, m N , respectively, e 



is the electron charge. In the limit Rn — Rs 3> £o, the 
superconducting pair potential can be approximated by a 
step function A(r) = A 9(i?;5 — r), where £o = hvp//S. 
is the coherence length, and vp is the Fermi velocity. 
Self-consistency of the pair potential is not taken into 
account, similar to the treatment given in 

Ref. Q At 

r = Rn, Dirichlct boundary conditions (i.e., an infinite 
potential barrier) are assumed, while at the NS interface 
r = i?5, the presence of a tunnel barrier is modeled by 
a delta potential V(r) = Uq 5(r — Rs). The energy lev- 
els of the system are the positive eigenvalues E of the 
BdG equation. In what follows, we consider the energy 
spectrum below the superconducting gap, < E < A . 
Rotational symmetry of the system implies a separation 
ansatz for the wave function as a product of radial and 
angular parts. We choose for the angular part the appro- 
priate angular-momentum eigenfunctions with quantum 
number m. Then the radial wave functions fm{r) satisfy 
a one-dimensional BdG eq. in the normal region: 
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with new dimensionless variables £ = r 2 /(2l 2 ) and e = 
E/(hui c ). The functions f^ir) are, respectively, the 
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electron and hole components of the radial Bogoliubov- 
de Gennes spinor. Here lo c — \eB\/m c s is the cy- 
clotron frequency, I — y/hc/\eB\ the magnetic length, 

v a = 2Ef ) /{hu c ), m± = &±rosgn(eB), £ s = fl|/(2Z 2 ), 
and sgn(x) denotes the sign function. After transforming 
the wave functions -> £ m ±/ 2 e~«/ 2 Eq. © 

results in a Kummer differential equation |l5j . and the 
ansatz for the wave function in the normal region {Rs < 
r < Rn) can finally be written as 
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where M{a,b,x) is Kummer's function [J^, and £y = 
R 2 N /{21 2 ). These wave functions satisfy the Dirichlct 
boundary conditions at r — R^, i.e., ifm^ {Rn) — 0, 
and the following symmetries hold: ifm' (r, e, £?) = 



(r, -£, --B) and </5m :±) (r, e, -B) 



(N,±), m 

<pL m '{r,s,B), 



where the dependencies on e and -B are emphasized for 
clarity. 

In the superconducting region r < Rs, the ansatz for 
BdG wave functions is given by |9j: 
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where 
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(r) = J m {q±r), q± = ky y/l± vrj, rj 
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Ao/(^TVAg-S 2 ), and J m (r) 
is a Bessel function of order to. These satisfy the sym- 



metries ipm \r, e) 



(S, 



and 7_ =7+- 
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The four coefficients a±,c± in Eqs. I|3a|l and JIJ are 
determined from matching conditions at the interface of 
the NS system [?| . These yield a secular equation for the 
eigenvalues e of the NS system for fixed mode index to. 
Using the fact that the wave functions tpm' + ^ given in 
Eq. (|3a|l are real functions and the symmetry relations 
between the electronic and hole-like component of the 
BdG eigenspinor, the secular equation can be reduced to 



Im{ 7+ D( n +)( e ,B)D( n -)( £ ,B)} =0, (5a) 
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and D { m \e,B) = 



D (+) (- 



e,-B) 



Here Z = 



(2ton/?i 2 ) Uq is the normalized barrier strength, and the 
prime denotes the derivative with respect to r. All func- 
tions are evaluated at r = Rs ■ The energy levels of the 
NS systems can be found by solving the secular equation 
(|5a|l for e at a given quantum number to. The secular 
equation derived above is exact in the sense that the usual 
Andreev approximation is not assumed [l6j |. An analo- 
gous result was found previously for zero magnetic 
field where the wave functions (pm are different. 

We now turn to the semiclassical treatment of the sys- 
tem. For simplicity, we assume a perfect NS interface, 



i.e., Z = 0, Ef ] = Ef' ms = m N . We follow the 
method developed in Ref. |9j, i.e., wave functions which 
appear in Eq. I|5b|l are approximated semiclassically. To 
construct these wave functions in the N region, one can 
use the standard WKB technique (see, e.g., Refs.H^ and 
ll^l for the radial Schrodinger equation with radial poten- 
tial given by V^t\r) = (t 2 /2 - m±) /r 2 . Here r = r/l. 
The four turning points (two each for the electron and 
the hole) can be obtained from Vm (r) = e. This yields 



p(N) 



T \2 = V 2 0± + m ±) T ^\Jv±{v± + 2to±) 
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where the sign in front of the second term under the 
square root distinguishes between the first and second 
turning points for both the electron and the hole, and 
v± = iso±2e. Note that r-j^ < r 2 , and the turning points 
are real if either m± > or v± > 2m± for m± < 0. For 
the electron and hole parts, the cycl otron radius g± and 
the guiding center c± are given |l9j by g± = I y/u ± 2e 
and c± = I yjv ± 2e + 2m±, respectively. 

The relative position of the turning points compared 
to ts = Rs/l an d tn = Rn/1 enables a classification of 
possible classical orbits which is summarized in Table [I] 
Orbits of type A\ correspond to the Landau states (or cy- 
clotron orbits) , while A2 are the so-called skipping orbits 
(or whispering-gallery modes discussed, e.g., in Ref. |^). 
In both cases, the orbits do not touch the superconductor 
and, hence, electron and hole states are not coupled. The 
other four types of orbits reach the NS interface. In case 
of type £?i , electron and hole alternately Andreev-reflect 
at the NS interface without ever touching the boundary 
of the N region. For type B2, the orbits reach both the 
inner and the outer circles delimiting the N region. Fi- 
nally, for types C\ and C2 , either the electron or the hole 
reaches the outer circle. 

In the S region, we approximate the wave function in 
the same way as in Ref.|9j. Substituting the corresponding 
WKB wave functions and their derivatives into the secu- 
lar equation l(5ajl and assuming Rs ^> £0, we obtain, after 
tedious but straightforward algebra, the following quan- 
tization condition for the semiclassically approximated 
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type of orbits 



Ai 



A 2 



Bi 



B, 



Ci 



Ci 



conditions 



TS < T 1 
T 2 < TJV 



TS < T x 
TJV < To 



T x < TS 
To < TAT 



Tj < T S 

t n < T 2 



T x < T S 
T2 < TJV 
TJV < T? 



Tj < T S 
T 2 + < TJV 
Tjy < T 2 ~ 



TABLE I: Classification of orbits. The solid/dashed lines correspond to the trajectory of an electron/hole. For the value of 
turning points Tf^, see Eq. @. 



energy levels: 



$m(ff) = n + fx 



(7) 



Here n is an integer, and the phase $ m (e) and the Maslov 
index [i are given in Table ILT1 The radial action (in 
units of h) of the electron and the hole between t\ and 
Ti reads 



5i ± )(T 2 ,T 1 ) = eL ±) (^^)-eL ±) (^r 1 ) , 



where 



-(v± 4- m±) arcsin 



v± + m± — t 2 /2 



\m± \ arcsin 



\Jv±(v ± + 2m±) / 

t 2 (v± + m±) — 2m± \ 
t 2 \Jv±{v± + 2m±) J 



(8b) 



Note that the value of the Maslov index comes out di- 
rectly from our semiclassical calculation. This result 
can be interpreted as follows. For orbits of type A\, 
the quantization condition can be simplified to = 
n + \ + \ {w>-± — |m±|), which coincides with the quan- 
tization of the electron/hole cyclotron states of a normal 
ring in a magnetic field. These are the familiar Landau 
levels. For orbits of type A2, the radial action between 
the boundaries of the classical allowed region (t^ and 
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Ai 


Sm ^ (t 2 , T-, ) 
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A 2 




3 
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Si 


Si+\T+,r s ) - sL\t-,t s ) - iarccos^ 





B 2 


Sm (t n , t s ) - S£> (t n ,ts) - \ arccos ^ 





Ci 


Sm } (t n ,t s ) - S ( m ] (T~ , ts) - - arccos ^ 


1 

4 


C 2 


Sin^T^jTs) - Sm (tjv,t s ) - i arccos 


1 

4 



TABLE II: Quantization conditions for the different orbits. 
See also the text. 



tjv) is equal to n + fi, where /i has a contribution j 
from the soft turning point (at r = 7^ ), and i from 
the hard turning point (at r = Tjy), resulting in an over- 
all /i = 3/4. For cases B\,B2,C\,C2, Andreev reflection 
takes place at the NS interface, resulting in an additional 
phase shift — - arccos^-. The action/Maslov index for 
the hole is —1 that of the action/Maslov index of the 
electron between the same boundaries. The conditions 
for the appropriate boundaries of these types of orbits 
can be obtained from Table |U The Maslov index is zero 
for cases B\ and B2 because the hole contribution cancels 
that of the electron. At the outer boundary for type C\ , 
there is a hard turning point for the electron and a soft 
turning point for the hole, resulting in /i = | — \ = \- 
Similarly, for C2, we have fi — — 5 + 1 = — \- 

In numerical calculations, it is convenient to use the 
following parameters that are suitable for characteriz- 
ing the experimental situation: k^Rs, where kp is the 
Fermi wave number, $ m i S s = BR 2 s ir is the missing flux 
due to the Meissner effect, Rs/Rn and A /Ep. Fig. Q] 
shows the comparison of the energy levels from the ex- 
act quantum calculation with the semiclassical results for 
two different systems. One can see that the semiclassical 
approximation is in excellent agreement with the exact 
quantum calculations. In Fig. the magnetic field is 
high enough for the appearance of Landau levels show- 
ing no dispersion as function of m. A small difference 
between quantum and semiclassical calculations occurs 
at the border of the region of the (cyclotron) orbits A\ . 

The experimental situation of Ref. 3J corresponds to 
the low-magnetic- field limit. There the cyclotron radius 
is large compared to Rn, and only orbits of type A 2 and 
B2 exist in the semiclassical approximation. Therefore, 
only these two types contribute to the free energy and, 
ultimately, to the susceptibility. In a simplified model, 
these orbits have been included in Bruder and Imry's 
theoretical study 0. Thermodynamical quantities such 
as the magnetic moment or the susceptibility can be de- 
termined from the energy levels of the system |2(J. How- 
ever, to fully explain the experimental results |3j, one 
needs to extend the work presented in this paper. For 
example, the Meissner effect can be included in a simi- 
lar way as in Ref. 0. The energy levels above the bulk 
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FIG. 1: Exact (crosses) and semiclassical (+ signs) energy levels (in units of Ao) obtained from Eqs. I5al and J7J as functions of 
quantum number m. In the left/right panel, k F Rs = 24.0, 18.0, "I> m i ss = 7.2, 4.05, and Rs/Rn = 0.29981, 0.14996, respectively. 
For both cases, Ao/Ep = 0.1. The solid lines represent the border of regions in the e,m plane where the different types of 
orbits arise. These lines can be obtained from conditions given in Table □ For easy reference, the different types of orbits are 
shown in the corresponding regions. 



superconducting gap {E > Ao) can be obtained by an- 
alytical continuation of the secular equation J5aj. One 
can expect a negligible effect from the roughness of the 
NS interface if the amplitude of the roughness is less than 
the wave length of the electrons 0, . 

In conclusion, we presented a systematic treatment of 
an experimentally relevant Andreev billiard in a magnetic 
field using the Bogoliubov-de Gennes formalism. An ex- 
act secular equation for Andreev-bound-state levels was 
derived that we evaluated both numerically and using 
semiclassical methods. In particular, a classification of 
possible classical electron and hole orbits in arbitrary 
magnetic fields was presented. This provides a useful 
starting point for successive studies of thermodynamic 
properties because it is possible to obtain the free energy 
of an NS hybrid system from the quasiparticle energy 
spectrum. Such an analysis may shed further light on the 
origin of the recently observed paramagnetic re-entrance 
effect and opens up a whole arena of new possibilities to 
study Andreev billiards in magnetic fields. 
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